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Abstract 

We construct a model of dark energy with a technically natural small contribution 
to cosmic acceleration, i.e. this contribution does not receive corrections from other 
scales in the theory. The proposed acceleration mechanism appears generically in the 
low-energy limit of gravity theories with violation of Lorentz invariance that contain a 
derivatively coupled scalar field Q. The latter may be the Goldstone field of a broken 
global symmetry. The model, that we call 0CDM, is a valid effective field theory up 
to a high cutoff just a few orders of magnitude below the Planck scale. Furthermore, 
it can be ultraviolet-completed in the context of Hofava gravity. We discuss the ob- 
servational predictions of the model. Even in the absence of a cosmological constant 
term, the expansion history of the Universe is essentially indistinguishable from that 
of ACDM. The difference between the two theories appears at the level of cosmologi- 
cal perturbations. We find that in 0CDM the matter power spectrum is enhanced at 
subhorizon scales compared to ACDM. This property can be used to discriminate the 
model from ACDM with current cosmological data. 
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1 Introduction 

Current cosmological observations indicate that the Universe recently entered a phase of 
accelerated expansion (see e.g. pLj and references therein). The driving force responsible 
for this phenomenon is normally called dark energy (DE) and its nature remains one of the 
major puzzles of contemporary physics. The standard explanation invokes a cosmological 
constant (CC), and the corresponding model - ACDM - manages to fit the data rather well 
[T] . However, on the theoretical side this explanation faces a serious problem. The measured 
value of the CC is many orders of magnitude smaller than the one expected on theoretical 
grounds. The latter is estimated as the value of quantum corrections to the vacuum energy 
plus the contribution related to possible phase transitions in the early Universe and is set 
by the energy scales of elementary particle interactions. This is the essence of the famous 
old Cosmological Constant Problem |2]. 

There exist two approaches to address this problem. The anthropic explanation assumes 
that the CC is indeed behind DE and claims that the observed value of the CC is small 
because the very existence of observers is possible only in a universe with small CC. The 
dynamical realization of this idea invokes the concept of a landscape of many vacua with 
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different values of the potential energy and transitions between them. Thus the underlying 
theory must be necessarily very complex, and it is a matter of ongoing debate whether this 
picture possesses any predictive power [3l HJ |5l [6l [7] . 

An alternative to the anthropic reasoning is to declare that DE is not related to the CC. 
For this mechanism to work, one first needs to get rid of the CC contribution to acceleration. 
One may speculate that there is some deep (presumably nonperturbative) mechanism that 
reduces the value of the CC, defined as the value of the total potential energy of all the 
fields at the absolute minimum, to zero [H |9l [10] . In this picture the observed accelerated 
expansion of the Universe must either be produced by some dynamical component of the 
matter sector or be due to modification of the laws of gravity at the cosmological distances 
leading to theories of dynamical DE (sometimes referred to as quintessence). Many models 
of this type have been constructed, see [11] for a review. However, a lot of them suffer from 
one or a few of common problems (which are often just other faces of the CC problem). The 
simplest models of quintessence involve very light scalar fields with the mass smaller than the 
present Hubble parameter. Generically, these masses are unstable with respect to quantum 
corrections which reintroduces a fine-tuning problem. With the additional assumption of an 
approximate global shift symmetry of the scalar field (in other words, assuming that this 
field is a pseudo-Goldstone boson) its mass can be protected from perturbative corrections 
[T2] . However, the VEV that breaks the symmetry corresponding to the pseudo-Goldstone 
field is required to have the value of order or above the Planck scale, which is believed to 
be problematic when non-perturbative quantum gravity effects are taken into account \13\ . 
Finally, from the purely observational perspective a drawback of this last class of theories 
is that, except for fine-tuned cases, their predictions are practically indistinguishable from 
those of ACDM which makes impossible the discrimination of these theories from the CC. 
From this perspective, the models of DE involving non-linear potentials and/or kinetic terms 
for the fields and their interaction with dark and ordinary matter are more interesting due to 
their non-trivial phenomenology. However, the non-linear structure of the potentials makes 
these theories strongly coupled at a low cutoff scale (typically of order 10 eV). Thus, these 
models can be understood only as effective theories below this cutoff, raising the issue of a 
missing ultraviolet (UV) completion at low energy scales. 

In this paper we present a simple model for dynamical DE that avoids the aforementioned 
problems. Assuming a vanishing CC, the smallness of the parameter setting the value of 
DE density in the model is protected by a discrete symmetry. The model is a valid effective 
field theory up to a very high scale just a few orders below the Planck mass. Remarkably, 
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above this scale the model admits a UV completion in the form of the Hofava's model |14j 
(more precisely, its consistent non-projectable extension fTEl [TBI [T7] ) which is a candidate 
for a renormalizable theory of quantum gravity. Finally, while the expansion history of 
the Universe in our model is basically indistinguishable from that of ACDM, the growth 
of cosmological perturbations is different, which can allow to discriminate this model from 
ACDM using the structure formation and weak lensing data0. 

The price to pay for these nice features is to admit violation of Lorentz symmetry in the 
gravity sector. This possibility is interesting on its own right and has received a great deal 
of attention recently, the prototypical models describing the possible patterns of Lorentz 
symmetry breaking being the Einstein-aether [19] and the ghost condensate [20] theories. 
A major issue in this kind of theories is the possible transmission of Lorentz breaking to 
the fields of the Standard Model (SM) which would be unacceptable as Lorentz invariance is 
tested to a very high precision in the SM sector [21j. Nevertheless, this problem is completely 
unrelated to DE and may be addressed using some dynamical mechanisms [221 |23l |23]. Once 
the emergence of Lorentz invariance at low energies for the SM is achieved, there is no extra 
requirement to make our proposal for DE technically natural. 

It is worth mentioning another interesting property of the model we are going to present 
that may shed some light on the old CC problem (vanishing of the vacuum energy). In the 
absence of a CC and matter (including dark matter) it admits two solutions, corresponding to 
Minkowski and de Sitter space-times. The latter solution is stable while the former is unstable 
with respect to long-wavelength perturbations. For that reason, the Minkowski solution is 
never reached during the cosmological expansion and instead the system approaches the de 
Sitter branch. In this sense the cosmology of our model is self-accelerating (cf. [25]). The 
instability of the Minkowski vacuum is a purely infrared phenomenon and is cut off at certain 
finite momentum. In other words, the Minkowski vacuum is perfectly well-defined from the 
point of view of the UV physics. This is consistent with the previous idea that the existence 
of such a vacuum may be required by some unknown principle which assures the cancellation 
of CC. In this paper we will not pursue the discussion about what kind of principle this could 
be. 

This paper is organized as follows. In Sec. |2]we present the model and describe the cosmo- 
logical solutions. In Sec. |3]we make a preliminary analysis of perturbations in the expanding 
and Minkowski backgrounds and show that in the latter case long-wavelength modes are 

non-trivial feature of the model is parametric enhancement of the signal in the structure formation 
that allows to have sizable effect even after satisfying stringent Solar system constraints. 
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unstable. In Sec. |4]we turn to the detailed analysis of the cosmological perturbations in the 
Friedmann - Robertson - Walker (FRW) universe and derive the set of linearized equations. 
These equations are solved analytically in Sec. [5] in various approximations. The reader who 
is interested in the final results can skip this section and go directly to Sec. Elwhere we report 
the results of the numerical integration of the linearized equations. Section [7] is devoted to 
conclusions. Appendix |A] contains some details about our numerical procedure. 



2 GCDM model for cosmic acceleration 

We start by considering theories where local Lorentz invariance is broken by the presence 
of a unit time-like vector u^. We will restrict to the case where the vector is orthogonal to 
a set of space-like 3-dimensional surfaces that foliate the space-time. Parameterizing these 
surfaces as the constant levels of a function f{x), can be written in the fornix, 

u,= , . (1) 

Note that is invariant under reparametrizations 

'P^<f = fif), (2) 

for an arbitrary monotonic function /. The most general action describing the coupling of 
Up to gravity and containing a total of two derivatives acting on Up has the fornjfl, 

r2 



M r 

S[EHx] = J d^x,/^{R + K^''^^VpU''V.u^) 



(3) 



where R is the space-time curvature, Mq is a mass parameter related to the Planck mass 
and 

= PS^pS: + \5i^S; + au^u'^g^p , (4) 

with a, /3, A free dimensionless parameters. This model has been proposecfl in [17] (see also 
[T8] ) and was dubbed 'khrono- metric' because it introduces the notion of a preferred global 
time coordinate set by the field ^p{x) - the 'khronon' field. As discussed in [17], Eq. ([3]) 
describes the low-energy limit of the consistent extension |T5l [16] of Hofava gravity [1^ 



^Our convention for the metric signature is (+, — , — , — ). 

■^We use subindices inside square brackets to distinguish the quantities referring to the different sectors 
of the model. 

''Note that the parameter A in corresponds to A' in the notations of [T7] . 
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The latter has been proposed as a power-counting renormahzable model of quantum gravity. 
In the case when the parameters a, /3, A are of the same order the action gives a valid 
effective field theory up to the scale Acutoff = min{Mo-y/a, Mq}. Note that when written in 
terms of the khronon field (p the action (l3l) contains four derivatives and thus gives rise to 
a fourth order equation of motion. However, the structure of the higher derivative terms is 
such that, for a proper choice of the time variable, only two of these derivatives are temporal 
[17] . Thus the khronon field contains only one propagating degree of freedom. In particular, 
no instabilities, that are in general associated with higher derivative Lagrangians, arise. 

Experimental data impose constraints on the parameters a, /3, A. For general order one 
ratios between these constants the strongest bounds come from the constraints on the pa- 
rameters af 012^^ of the parametrized post-Newtonian (PPN) formalisnj^ and require 
a, /3, A to be smaller than 10~^ [17j. However, for the special case a = 2(3 both af'^'^, af^^ 
vanish and one is left with a much weaker bound a, /3, A < 0.1 coming from the dynamics of 
the Big Bang Nucleosynthesis (BBN) p7j and the emission of gravitational waves in binary 
systems |28j. Another special case is /3 = 0, a = A. Then the PPN parameter a2^^ vanishes, 
while the constraint on af gives a, A < 10""^. 

If we do not insist on the existence of a known UV-completion and are content with 
an effective theory with high cutoff, we can relax the condition ([T]) and consider a general 
unit time-like vector. This would lead us to the Einstein-aether theory [19] (see |29] for a 
review). The action for the general aether has the same form as ([3]) with an extra ternj^. 
Physically, the Einstein-aether model differs from the khrono-metric theory by the presence 
of transverse vector modes among the perturbations. However, these modes do not affect 
the homogeneous cosmological expansion, nor the linear evolution of scalar cosmological 



perturbation^!!. Thus, the analysis of the present paper applies without changes to the 
general aether model. Just for clarity we concentrate on the hypersurface-orthogonal case 
corresponding to the khrono-metric theory. 

To the system ([3]) we add a scalar field G with the symmetry under the shifts 

e + canst . (5) 

The precise origin of this field is not important for our purposes. However, an example that 
is useful to keep in mind is a Goldstone boson corresponding to a spontaneously broken 



^ These PPN parameters describe effects of local Lorentz violation and are zero in general relativity. The 
experimental bounds are [25]: |af^^| < 10~^, |q;|'^^| < 10~^ . 

^For a hypersurface-orthogonal aether this term can be written as a combination of those present in ([3|) . 
^The transverse modes can lead to interesting effects in the CMB polarization [50] . 



6 



global symmetr > ^. For this reason, we will refer to the field G as "Goldstone field" in the 
rest of the paper. In contrast to the pseudo-Goldstone quintessence models [12], we assume 
the symmetry ([5]) to be exact. The interaction between the Goldstone and the khronon must 
be invariant under ([2]) and ([S]). Keeping only the operators that have dimensions up to 4 
and thus are relevant for physics at long distances, we obtain the following general actioij^ 
for 9: 

^[0] = / d x^-g ( ^ + ^ + u^Of,Q \ . (6) 

Here the coupling constant x is dimensionless, while /i has dimension of mass. It is important 
to notice that the /i-term is the only one that is not invariant under i— )■ —0. Thus, assuming 
that the UV completion of the theory respects this discrete symmetry, the value of /i is stable 
under radiative corrections. 

The physical effect of the two khronon-Goldstone couplings in ([6]) is quite different. Let 
us consider the dynamics of the Goldstone field in the background with fiat metric and 

= (1,0,0,0). Then the third term in ([6]) becomes a total derivative and does not affect 
the dynamics, while the second term modifies the propagation velocity cq of the Goldstone 
field: 

" 1 + X • 

We will see in a moment that the effects of the khronon-Goldstone coupling described by 
the third term become essential in curved backgrounds and modify the dynamics of both 
fields at cosmological scales. A coupling of this form was introduced in fSlj in the context 
of infiation in the Einstein-aether model. A similar coupling was also considered in [32]. 
We will study the cosmological evolution in the model with the action 

S = S[EHx] + ^[B] + S[m] , (7) 



^Note that the scale of this putative symmetry breaking wih be completely irrelevant for the analysis in 
this paper (provided that it is higher than the temperature of the universe at the epochs of interest, which 
is roughly a few eV). In particular, it can be safely below the Planck mass. 

^ One could consider including in the action operators 

that formally have dimension 3. We omit them for two reasons. First, these operators can be forbidden 
by imposing the symmetry under simultaneous change of the signs of Q and u^: i— > — , i— )• — u^. 
More importantly, these operators actually have dimension higher than 4 when expressed in terms of the 
canonically normalized khronon perturbations, and thus do not affect the long distance physics. 
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where S[ehx] ^[e] have been introduced above and S[m] stands for the action of matter, 
which includes baryonic matter, dark matter and radiation. We assume that S[m] has the 
standard form with matter interacting minimally with the metric Qf^u- In particular, we 
assume that there are no direct couplings between the matter sector and the khronon or 
Goldstone field£°|. We will refer to this model as 6CDM. 

Introducing the Ansatz for the most general spatially homogeneous solution 

ds^ = N'^{t)dt^ - a'^{t)6ijdx'dx^ , (8a) 

^ = ^it), e = e{t), (8b) 

into the actions ([3]), ([6]) we obtain 

r / ^3 



M2(6 + 3/3 + 9A) r 4 d^a f f a • 2 , 2 3A 
S[EHx] + S[e] = J d^x— + J d^xi^-^^e^ + fx^a^Q 

Note that the field if has dropped out of the action, which is a consequence of the invariance 
under ([2]). Varying the action with respect to 6 and we obtain the cosmological equations. 



^(j^r^ , ,.2.3 

dt\Ncl 



e + /i'a^ = , (9) 



^2 



where H = d/{aN) is the Hubble rate and 

G 1 

^cosm 8^^2(1 + + 3A/2) ■ 

On the r.h.s. of the Friedmann equation (fTOj) we added the contribution pm of ordinary 
matter. Note that the expression for the gravitational constant Gcosm entering the Friedmann 
equation is corrected by the presence of the khronon compared to the standard GR formula 



^''Tliis must be considered as a simplifying assumption. For the ordinary matter the direct couplings to the 
khronon field would imply violation of the Lorentz symmetry within the SM sector and is tightly constrained 
by the experimental data. Once these constraints are satisfied, the effects in cosmology are negligible. On 
the other hand, for dark matter we do not have such strong observational evidence of Lorentz invariance. 
Allowing for direct coupling between dark matter and the khronon may lead to interesting effects that will 
be reported elsewhere [33]. As for the coupling of SM / dark matter to the Goldstone, due to the symmetry 
([5]) it must involve derivatives. Apart from possible kinetic mixing with a scalar gauge singlet that may be 
part of dark matter, this would produce higher order operators irrelevant for long-distance physics. 
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G = (SvtMq) ^. As discussed in [l5l[T7], this differs from tlie expression for tlie gravitational 
coupling 

r 1 

^ 87rM2(l-a/2)' 

appearing in the Newtonian potential of a localized source. The bound [27] \Gcosm/GN~i\ < 
0.13 from BBN thus requires the parameters a, /3, A to be smaller than 0.1. 
The equation ([9]) is easily solved with the result 

(11) 

where C is an integration constant. Due to the cosmological expansion the ratio Q/N is 
attracted to the constant non-zero value —(i'^Cq. On this attractor solution the field G 
produces the same contribution into the Friedmann equation (ITOl) as the vacuum energy 
density 

Vo,eff = . (12) 

In this way the model leads to accelerated cosmology even in the absence of the true cos- 
mological constant. Due to the stability of the parameter /i under radiative corrections it is 
technically natural for the effective vacuum energy density ( fT2l) to be small. 

Notice that the general solution ( fTTl) . once introduced in the Friedmann equation ( JTOj) . 
also produces a contribution scaling like a~^, that one may be tempted to interpret as dark 
matter. However, this contribution is always accompanied by the stiff term proportional 
to C^a~^ that dominates over the dark matter-type contribution. To comply with the ex- 
perimental constraints the stiff term must be small. We will set C = in what follows. 
Thus the 0-field cannot be responsible for dark matter, which one has to add as a separate 
component. This will be included in the part of the theory described by S^m] in ©• 

It is worth stressing that we do not present any mechanism that could lead to vanishing 
of the vacuum energy and thus do not pretend to solve the old CC problem. However, we 
point out the following interesting property of the BCDM model that may shed some light 
on the resolution of this problem. In the absence of matter, pm = 0, the system of equations 
(E]), (fTOl) possesses a static solution with fiat metric, 

N = a = const, 9 = 0. (13) 

However, this solution is never achieved if the evolution starts from a configuration with 
Pm 7^ 0. Moreover, in the next section we will see that the Minkowski solution is unstable 



9 



with respect to long-wavelength perturbations. On the other hand, at small scales the 
instability is cut off and from the viewpoint of the short distance physics (IT^ is a valid 
background. Thus one may hope to find a mechanism that would neutralize the CC in 
this background and ensure the existence of a (unstable) Minkowski vacuum. However, the 
dynamics drives the universe not to Minkowski but to the de Sitter attractor. 



3 Preview of linear perturbations 

As we saw in the previous section, at the level of homogeneous cosmology the attractor 
solution in the BCDM is indistinguishable from ACDM. Then it is natural to wonder if the 
model leaves any distinctive signatures at the next level of approximation, namely in the 
evolution of linear cosmological perturbations. The analysis of the linearized perturbations 
is also indispensable to verify the stability of the accelerated solution. 

Before going to the complete analysis which we postpone to the subsequent sections, let 
us consider a simplified setup where one neglects back-reaction of the khronon and Goldstone 
perturbations on the geometry. Formally this corresponds to taking the limit Mq — ?■ oo. As 
will be clear later, we must at the same time keep the scale 



(14) 



fixed. This implies that in the limit of interest a — > 0. We also assume that the ratios P/a 
and A/a are fixed and of order one in this limit. 

We start with the analysis of perturbations in Minkowski space-time. One writes 



V9(t,x) = Lp{t) + X(i,x) 

e(t,x) = e(t) + e(t,x) 



(15a) 
(15b) 



where ip{t), G(t) are the background values of the fields and ^ are the perturbations. Note 
that due to the reparametrization symmetry ([2]) we can always set 



<^(t) = t . 



(16) 



Also for the Minkowski solution 0(t) = const. Expanding the actions 
order in fluctuations we obtain 



[S[EHx] + 'S'[0]]^^^ ^ J ^ 



X 



to quadratic 







(17) 
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where A = didi is the spatial Laplacian, and we have introduced the notation 



c = . 18 

^ a 



Throughout the paper we will assume that c^, as well of order one. The action 

ffT7|) yields the equations of motion, 



= 0, (19a) 

4- - Ae - 12' Ax = . (19b) 

The overall Laplacian in the first equation can be cancelled out for any inhomogeneous 
configuration. Performing the Fourier decomposition 

X, e oc 6'^'^*+''^'' (20) 

we find the dispersion relations of the propagating modes, 

1 



(4 + 4,)k^±^Hci-cirk^ + ^k^- 



a 



(21) 

This expression simplifies in two regimes. At k larger than the critical momentum 

h = /iV^a , (22) 

the fields x C decouple and describe two modes with linear dispersion relations 

col = c^e , (23a) 
Lol = c|fc2 . (23b) 

Note that both modes are stable. In the opposite case k Eq. (^1^ takes the form, 

cul = ±cekck . (24) 

This dispersion relation is non-analytic: the frequency is proportional to the square root 
of the momentum. Clearly, one of the modes has purely imaginary frequency and thus 
exponentially grows with time. This signals instability of the Minkowski background with 
respect to long-wavelength perturbations. As we are going to see, this instability disappears 
in the expanding universe. 
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As the next exercise we consider the evolution of the X'l ^-perturbations in an external 
FRW metric. We will work in the conformal time which corresponds to setting = a in 
f l8a|) . and consider a general time-dependence of the scale factor. We expand again as in 
f lT5]) where we fix f[T^ . The background configuration G)(t) satisfies Eq. flTTl) . To simplify 
the calculations we concentrate on the attractor solution 

9 = -a/i^cl . (25) 

After a straightforward computation we find the quadratic action 

[S[EHx] + 'S'[e]]^^^ = j d^xa^ 



{Mln{l -B) + MX(1 + B) + aV4)^^ (26) 



\2 

,2^2 



+ 771" n ^/^ CediXC'i^ 



where we introduced the notations 

fl = ^, (27) 

a 

H = - = Ha, (28) 
a 

Note the appearance of a term with two space derivatives of x which was absent in the case 
of Minkowski background. The action ( l26l) yields the following equations, 

X + 2nx- clAx + [nil -B)+ n\l + B) + a'c%kl]x + = , (29a) 

e + 2ni - c| - a/i^cl Ax = , (29b) 

where in the first equation we have cancelled the overall Laplacian. The crucial difference 
from Eqs. f[T9|) in the Minkowski case is the presence of an effective mass term for the field 

™?„^^Mi,^^^ + 4^?. (30) 

This expression deserves two comments. First, during the epoch of primordial inflation one 
can neglect the last term in (!30|) and also use T-i ^ "H^. This gives ^e// ~ '^■^infi 
the effective khronon mass is of order the Hubble parameter at inflation and coincides with 
the effective mass of a conformally coupled scalar field. As a consequence, the khronon 
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perturbations are not generated during inflation. This agrees with the result of [34j that no 

isocurvature perturbations associated with the khronon are produced during inflation. 

Second, for certain choices of parameters the effective mass square f pOj) can become 
negative during the radiation and matter dominated epochs. This will lead to the growth of 
the khronon perturbations with wavelengths larger than the horizon size at the corresponding 
epoch. It is worth stressing that this growth is under control, its rate being proportional to 
the Hubble parameter. Even if this regime is definitely interesting from the phenomenological 
point of view, it corresponds to a rather special corner of the parameter space. In the present 
work we focus on the generic stable case. 

Let us now analyze the behavior of subhorizon modes with large momenta and frequen- 
cies, k, uj ^ H. In this case we can neglect all 7^-dependent terms in Eqs. (129|) and perform 
the Fourier decomposition f l20|) . However, the presence of a non-trivial background 6 still 
modifies the equations for the perturbations, eliminatin g a ny instability for the momenta 
inside the horizon. Explicitly, the dispersion relations areo 



For large momenta k ^ kc one recovers the linear dispersion relations ( l23l) . While at /c ^ fcc 
we have 



We see that, in contrast to the case of Minkowski, both modes are stable (cf. (12^ ). One 
of them possesses a frequency gap, while the other has a dispersion relation with quadratic 
dependence of the frequency on momentum. This implies that the latter mode has low 
propagation velocity c_ ~ k/k^ ^ 1 and one expects it to cluster in the gravitational 
potential wells enhancing the growth of the cosmological perturbations. We are going to see 
below that this expectation is correct. 

It is worth stressing that the analysis of this, as well as the subsequent sections, applies 
without change to the longitudinal sector of linear perturbations in the case when the khrono- 
metric model is substituted by the Einstein-aether theory to describe the Lorentz breaking. 

"'^"'^Here w and k should be understood as physical frequency and momentum which corresponds to putting 
a = 1 in the formulas. 




(4 + ci)e + c%ki ± ^{{ci + ci)k^ + c%kiY-Acicik^ 




(31) 



(32) 
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4 Linearized equations in FRW 



We now turn to the derivation of linearized equations of motion for ©CDM including gravi- 
tational perturbations and matter fields. We start with the Einstein's equations 

-^M^^ - 2^'"'''^ ^ + ^[©Im'^ + ^M/ii.) , (33) 

where T^^q^^^, and T^m]^y are respectively the energy-momentum tensors of the khronon, 

Goldstonc field and matter. They are defined as the variation of the corresponding action 
with respect to the metric, 

2 5S^^ 

After a long but straightforward computation we obtain 

1 (34) 
+Vx{K(^^u^)) +a[axa^u^Ui, - 2a^V ^^uxu^) + a^,ay\ + -K^^^xu^g^u , 

where 

= u^VxUfj, , (35a) 
= K^^^xu" , (35b) 

and the round brackets denote symmetrization of indices: 

Similarly, for T^Qj^^,^ we find, 

T[e]M. = d^Qd.Q + 2{>€u''dxQ + /x') ^(A)© " u^u,{>i{u^dxef + fi\^dxe) 

'V^BSaO X, x^ 2 ^^^^ 

— + ^(«^A©)' + /xV^aO j^^, . 

Finally, for the matter sector we will include only two decoupled components correspond- 
ing to cold matter (cm) and radiation (7) in the perfect fluid approximation. Reducing the 
matter sector to this form is, of course, a crude simpliflcation but it is convenient for illus- 
trating the difference between 0CDM and ACDM. In this respect, whenever we refer to the 
ACDM model in the future, we will have in mind the case of GR with a CC and the previ- 
ous matter content. We leave more accurate analysis of the cosmological perturbations for 
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future research. Thus, for the energy-momentum tensor we take the standard hydrodynamic 
expression for each component 

T[m]^lu = iP[j] + Ph])^h]t^^h]iy —P[-t]gtJ.y + (P[cm] + P[cm\)'V[cm]^lV[cm]v , (37) 

where p[a], p[a] are the energy density and pressure of the different matter components and 
v\^a]n are their 4-velocities. 

To form a closed system the Einstein's equations must be supplemented by the equations 
of motion for the khronon, Goldstone field and matter. The first two of these equations are 
obtained by varying the sum of the actions ^ and ([HD with respect to ip and 6. This yields, 



pp^l 



(38a) 
(38b) 



where 



The equations of motion for matter follow from the conservation of the energy-momentum 
tensor of each component (recall that we assume the two components to be decoupled). 



V^Tr 



. 



(39) 



In fact, only one of the equations in (139|) is an independent equation, as the conservation of 
the total matter energy- momentum tensor follows from the Einstein's equations and fl38|) . 

We now expand the above expressions to linear order in perturbations around a FRW 
background. The dynamics differ from the standard case only in the scalar sector of the 
perturbations and we concentrate on this in what follows. For the khronon and Goldstone 
fields we take the representation ( fTSl) with the background values satisfying ( fT6l) . ( !25l) . The 
metric perturbations are chosen in the conformal Newton's gauge 

ds^ = a\t) [{1 + 2(j))dt^ - (1 - 2^)6ijdx'dx^] . 

Finally, for the scalar sector of the matter perturbations we write 

P[a] = P[a]it) + Spia] , P[a] = P[a]{t) + Sp[a] , f[a]j = djVla] . 
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Substituting these expressions into Eqs. ( 133|) . ( 136|) . (!37|) we obtain the expressions for 
the hnearized Einstein's equations in 0CDM, 

2A^ - 3'H{2 + aB)i) - aA^ - '^iM±p!I^^ + ^tAx + a(l - 5)?/Ax 

Mq 

u^a- aHSp[^] + 5picm]) _ , . 

^mI^ Ml - " ' ^^"^^ 

(2 + aB){i, + net.) + (/3 + A) Ax - ^^^^^^p^^H - = , (40b) 

3(2 + aB) + H{4) + 2^) + {2H + 7/^)0] + 2A0 - 2AV^ 

+a5(Ax + 2?^Ax)-^^=0, (40c) 

0-^-/3(x + 2?{x)=O, (40d) 

where B and "H are defined in (!27|) . (!28|) . These equations are respectively the (00), (Oi), the 
trace and trace-free parts of the {ij) Einstein's equations. In deriving fl40cl) we have used 
that the background metric satisfies the equation 

^ ' Ml Ml 

which can be obtained by taking the time derivative of ( [T0|) . Note that according to Eq. (]40dp 
for non-zero /3 the khronon field induces anisotropic stress and thus the gravitational poten- 
tials and are in general different. 

The equations of motion (138|) for the khronon and Goldstone field read at the linear order: 

x^mx- c\Ax + (^(1 - 5) + n\\ + B)Wc%kt)x^ 

= + 7/(1 + 5)0 + 5^, (41a) 
I + 2H^ - 4 A^ - afi^c^Ax = -a/i^40 - 3afihln(f) , (41b) 

where we remind the definitions ( ITSl) . (!22|) of C;^ and kc- Note that these equations differ 
from Eqs. ( |29l) valid in an external FRW universe by the appearance of the source term 
proportional to the metric perturbations on the r.h.s. 

The perturbed matter equations coming from (!39l) have the standard form: 

5p[,] + 3H(5p[,] + 6p[a]) - (P[a] + P[a]) + 3^^ = , (42a) 

V[a] + mv^a] + ———m - — — — - a0 = . (42b) 

P[a] + P[a] P[a] + P[a] 
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In the next sections we analyze the system of hnear equations (140|) . ( 14T|) . fH2|) . 



5 Cosmological perturbations: qualitative analysis 

We begin with a quahtative discussion of the evolution of cosmological perturbations in the 
0CDM model at radiation and matter domination. From Eqs. f l40l) it is easy to see that 
at these stages the modifications of the linearized equations of GCDM, as compared to the 
ACDM case, are proportional to the dimensionless parameters a, f3, A. As discussed in Sec. [21 
consistency with the experimental data requires these parameters to be small, A ^ 1. 
This allows to solve the system ( l40l) . ( 1411) . f l42l) with the following perturbative scheme. One 
first solves Eqs. (HOj) . (H2|) in the zeroth-order approximation neglecting all the contributions 
proportional to a, (3, A. Assuming adiabatic initial conditions, the solution is given by the 
standard adiabatic mode. For clarity we concentrate on the modes entering inside the horizon 
at the matter-dominated epoch. In this case the adiabatic mode is particularly simple: the 
two gravitational potentials (p and ip are equal and constant during radiation- and matter- 
dominated epochs, with a jump at the radiation-matter equality. 



At the second step of the perturbative scheme these expressions are substituted into Eqs. ( HTl) 
as the source for the fields x ^i-nd ^. Finally, the solution for x and ^ is inserted back into 
Eqs. (I40l) to produce the corrections 4'^-^'^ aiid (5 p[a] / pia]Y^^ to the observable quantities. 
Thus our task at the moment is to solve Eqs. (14T1) with the sources ( 143|) . 

A complete analysis should also include the evolution at the epoch of DE-domination. 
However, as this epoch started only recently it does not quahtatively affect the results. At 
the same time, all the following analytic estimates should be understood as giving only 
the qualitative picture. For example, they do not capture the transient behavior at the 
boundaries of the different dynamical regimes, which may be important due to finite duration 
of the corresponding epochs. We will neglect these effects in the rest of this section, and 
postpone their proper treatment to the numerical analysis (see Sec. [6]). 




(f).y = const , radiation domination; 
(pcm = , matter domination. 



(43) 
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k 



kc 



ki/2 

h/3 
Ho 

teq to t 

Figure 1: Time dependence of different physical scales relevant for the dynamics of cosmolog- 
ical perturbations: Ti (blue solid line), akc (magenta dashed line), y/l-iakc (green dash-dotted 
line). The time of the triple intersection of these curves is denoted by tx- The evolution of 
perturbations proceeds differently in the regions (a'), (a), (6), (c), {d). The time of radiation- 
matter equality is denoted by teq] to is the present time and Hq is the present Hubble rate. 
The present scale factor a (to) is set to 1. The beginning of current accelerated expansion 
roughly corresponds to the point where l-L starts growing. The definition of fcc, ^1/2, ^1/3 is 
given in the text. 

5.1 Different regimes of evolution 

One expects the behavior of a mode with a given comoving momentum k to be different 
depending on the relation betweeiJ^ k and various scales that can be constructed out of 
the coefficients in Eqs. (j^T]) . Clearly, two important scales are the Hubble rate l-L and the 

""^^In principle, the comparison of the different scales should include a dependence on the propagation 
velocities ce and c^. For simplicity we will consider the case where both quantities are close to one. 
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comoving value of the critical momentum akc- Below we will encounter a third relevant 
scale, the geometric mean of the first two, y/Tiakc- The meaning of this latter scale can be 
understood as follows. For modes with k = y/Hak^ the frequency (|32|) of the low-frequency 
branch is comparable to the expansion rate of the universe and therefore the modes with 
this and smaller k are significantly affected by the Hubble friction. The dependence of the 
three scales on time is illustrated in Fig. [1] One distinguishes five different regimes of the 
mode evolution: 

(a') superhorizon modes, k < Ti, radiation dominated universe; 

(a) superhorizon modes, k kT-L, matter dominated universe; 

(6) subhorizon modes with "H, akc < k; 

(c) subhorizon modes with "H, yjl-iakc < k < akc] 

{d) subhorizon modes with Ti < k < \/7lak^. 
For the following discussion it is convenient to normalize k to be equal to the physical 
wave-number of the mode at present by setting the today scale factor to unity, a{to) = 1. It 
is also useful to express k^ in terms of the present Hubble rate Hq. The latter is given by 
the expression 



where Qde ~ 0.75 is the dark energy density fraction. Recalling the definition (B^ we obtain 



We see from Fig. [T] that depending on k the different modes go through different regimes 
of evolution. Before studying these regimes case by case, let us briefly outline their broad 
features. 

The modes with k > kc go only through the regimes (a'), (a) and (b). According to the 
analysis of Sec. [3], in these regimes mixing between the x- ^^id ,^-perturbations is weak and 
both modes have linear dispersion relations with order one velocities. One expects that the 
effect of these modes on the growth of cosmological perturbations is at most of order a, (3, A, 
as it is the case for pure Einstein-aether and khronon theories [55] . 

The dynamics of the modes with ki/2 < k < kc are more interesting. Here 



These modes spend some time in the regime (c). We are going to see that in this range 
there are non-trivial effects due to x - ^ mixing and the backreaction onto the gravitational 





(44) 




(45) 
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potentials is enhanced by a factor a^^/^ compared to the short wave-length regime. Note 
that for a not extremely small, a > 10^^, the momentum ki/2 lies below the Hubble rate 
at the radiation-matter equality!^ Heq ~ IOOHq. Thus the modes with k ~ ki/2 enter the 
horizon at the matter-dominated epoch. 

The modes with ki/^ < k < ki/2 pass through all five regimes {a')-{d). The momentum 
fci/3 corresponds to the intersection of the lines T-L(t) and a{t)kc, see Fig. [H To determine 
it we note that for the interesting values of the parameters the intersection happens at the 
epoch of matter domination. During this epoch the scale factor is given by 

Q = ^cmt , ^cm = ^ ) (46) 

where f2cm ~ 0.25 is the matter density fraction. This gives for the time of the intersection 

t -( ^ \" 
y^cmH^kc J 

and hence, 

For these modes we will also find a relative enhancement of the backreaction onto the grav- 
itational potentials of order a"^/^. 

Finally, the modes with k < ki/^ go directly from the regime (a) to the regime {d) 
and do not have enough time to develop. Thus the backreaction of these modes on the 
gravitational potentials is expected to be 0(a). Note that for these modes to be observable, 
their wave-length must be inside the present horizon size, k > Hq. 

We now confirm the statements made above by the detailed analysis of the listed regimes. 

Regime (a'): k during radiation domination. 

As we are considering superhorizon modes, we can neglect all terms with spatial Laplacians 

in Eqs. (14T]) . Some care is needed with the last term on the l.h.s. of ( ]41bp as it is the only 

term in this equation containing the field we will check explicitly below that this term is 

indeed small on the solutions. 

^■^This relation is obtained as follows. At the matter-dominated epoch a cx and H = 2/t. Hence, 
Heq = Ho{to/teq) — Hoy/z^q + 1, wherc Zeq ^ 10** is the corresponding red-shift. 
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At the radiation-dominated stage the scale factor depends hnearly on time, 

a = A^t , A^ = y^Ho , (48) 

where fl^ ~ 10~^ is the radiation density fraction today. Substituting this and ( H3|l into 
Eqs. ( HTl) we obtain 



2 



A^tclkf.^ 1 + B 



X+ " 2 '^ e = ^— (49a) 
e + ^^ = -3AA|0^. (49b) 



This system has the solution: 

^ ^ 07^ ^ ^ ^7/i^c|0^t^ ^^^^ 

It is straightforward to check that on this solution the ratio of the term Ax in (]41bp to the 
term ^ is of order c%,k'^t^ which is indeed small for superhorizon modes. 

Note that despite the apparent growth of the fields x, i in (l50l) . the physical perturbations 
remain subdominant with respect to the background. Indeed, combining (150|) with the 
expressions (ITBi) . for the background we obtain the total values of the fields Lp and 6: 

<^ = t (1 + 0^/2) , 



,,2 2 /I f2 

We observe that the relative corrections are small. The solution fl50|) is equivalent to the 
time shift 

th^t (1 + 0^2) , 

that corresponds precisely to the adiabatic mode at radiation domination. 

Let us check whether the solution fl50|) is an attractor. To this end consider the solutions 
of the corresponding homogeneous equations (i.e. Eqs. (jlH]) with zero r.h.s.). Those are 
particularly simple under the assumption that the second term in the square brackets in 
f l49ap can be neglected. This is equivalent to the assumption 

l-L->akc, (51) 

which is always satisfied during radiation domination for reasonable choices of parameters. 
Then the solution of the homogeneous equations is 

^ 2{6 + B)fi^^^ 2(3 + 5)/x2 ^ + ' ^ ^ 

^hom = ^0 + y , (52b) 
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where 

-1 ± VI -8fi 

= 2 

and C*) -D_ are integration constants. 

Let us first discuss the contributions proportional to and C. As long as we concentrate 
on the field ^, the role of these contributions decreases with time compared to the adiabatic 
mode (l50|) . Moreover, the constant mode is irrelevant when we are interested in the back- 
reaction of ^ on the gravitational potentials. Indeed, the only place where ^ appears in the 
Einstein's equations is Eq. (140 ap and there it enters with a time derivative. On the other 
hand, the terms proportional to and C in the homogeneous solution f l52ap for x grow with 
time faster than the adiabatic mode and may become important at late times. Whether this 
happens or not depends on the initial conditions for ^ which set the values of .^o and C. A 
proper determination of these initial conditions would require analysis of the evolution of 
the system at the epoch prior to radiation-domination (i.e. at infiation and reheating) and 
also a precise model for the origin of the field G (e.g. if it exists during all the history of 
the universe, or appears as the result of a phase transition at some particular epoch). These 
issues are outside the scope of the present article. Nevertheless, let us give an argument 
showing that ,^o- and C-contributions into x are likely to be small. For concreteness we 
concentrate on the ^o-contribution. Consider the ratio of this contribution to the adiabatic 
mode ( 150|) . Omitting factors of order one and using ( 1T4|) . ( HHj) it can be written in the form, 

where Hinf is the Hubble rate at infiation and e is the infiationary slow-roll parameter. In 
passing to the second expression we have used the standard formula 



Gcosm Tj- 

e 

for the amplitude of the primordial perturbations in the potentially driven slow-roll infiation. 
The first factor in f l53|) is small due to the inequality (|5T]) . Let us estimate the second factor. 
Assuming that the fiuctuations of the ^-field are generated at the infiationary epoch, we have 
the estimate ^ Hmf. We will be interested in the values of a in the range from 10~^ to 
10~^. Finally, recalling that the measurements of the tilt of the primordial spectrum give 
an upper bound e < 10^^, we obtain that the second factor in ( l53i) can be at most < 10. 
Overall, we conclude that the product (153|1 is small in general. A similar reasoning applies 
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to the contribution in fl52ap that is proportional to C. We will neglect these contributions 
in what follows. 

Consider now the last two terms in f l52ap . The term proportional to always decays 
with time and can be safely neglected. However, the term can either decay or grow 
depending on the sign of the parameter B. Moreover, for B < —1 this term grows faster 
than the adiabatic contribution fl50|) . We have already encountered this growing regime in 
Sec. [31 where it was related to the possibility of having negative khronon mass squared. 
Below we will only consider the case B > —1 where the leading solution is given by the 
adiabatic mode fl50|) . 

Regime (a): k Kl-i during matter domination. 

During matter domination the scale factor is given by the formula (H6i) . Substituting this 
into Eqs. (HTl) and neglecting the terms with spatial Laplacians we obtain 



4 

X + jX 



4^ -.2 2 



cm ^0^c 



X + ^=^?^^^.™. (54a) 

yU^ t 



^ + 7^ = -6Acmt/i'40em • (54b) 



The adiabatic mode is also a solution and has the form 



It corresponds to the time shift 



Similarly to the previous case one can check that the term proportional to Ax on (141bp can 



be safely ignored. Also, as in the previous case, one can argue that the solutions of the 
homogeneous equations are irrelevant, apart from one, possibly growing, mode. This mode 
appears in the regime (15T|) and has the form 



-3 + VI - 245 
Xhom oct+ , r+ = . 

It grows faster than the adiabatic mode if B < —1. This growth is cut off once the mode 
enters inside the horizon or the inequality flSTj) is violated. Below we restrict to the stable 
case B > —1. 
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Regime (6): subhorizon modes with ak^ < k. 

We now turn to the evolution of the subhorizon perturbations. We shall concentrate on 
the modes that enter inside the horizon at the matter-dominated stage. Neglecting the H- 
dependent terms in brackets in Eq. (14 lap and substituting (H3l) on the r.h.s. we obtain, 



X + JX+ [elk' + A^t^kl] X + = ^^^<^cm , (56a) 

+ clk'^i + A,^e^i'c%k\ = -6A,^tfi'cl(f),^ . (56b) 

Deep inside the regime under study, k ^ akc, the previous equations have the approximate 
solution: 

^ = ?(i±^, ^ = -2(3 + Ha(i + B))^Wfc^. (57) 



c'ik'^t ' \ cl ) k"^ 

A, ^ /C 



This solution decays with time (in the case of the .^-perturbation it must be compared with 
the behavior G oc t'^ of the background). On top of fl57|) the fields x? ^ exhibit oscillations 
with frequencies u-^, (Eqs. ( l23l) ) and amplitudes decaying as These oscillations 

represent small corrections to the approximate solution (|57l) which disappear at large times. 



Regime (c): subhorizon modes with ^/l-iakc < k < akc- 

According to the results of Sec. [3l in this regime one of the branches of x^C perturbations 
has a frequency gap of order kc, see Eq. ( 13T|) . As long as we consider frequencies smaller 
than the gap, this mode can be integrated out. To obtain the equations for the remaining 
mode we use the following trick. First, we eliminate the explicit dependence of the equations 
on the parameter fi by introducing 

The second step is to take the limit 

fee — )■ oo . (58a) 

The physical meaning of this limit is that we focus on the modes with frequencies u much 
smaller than the gap (cf. fl32|) ). At the same time we must be careful to keep the frequency 
of these modes non-zero. From Eq. ( 132|) we read out a;_ oc k'^/k^. Thus to keep fixed 
f l58ap must be accompanied by 



A; — )■ cxD , k (X \/kc ■ (58b) 
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The rest of the quantities in the equations must remain finitcj. In this way Eq. fl56ap reduces 
to 



X + e = 0, (59a) 



I A +2 2 
A +2„2 1,2 ^ ^^"^'^ 



where we have kept the main subleading contribution represented by the first term in the 
brackets. Substituting this into fl56b|) we obtain, 



•■ 4 ~ c^k'^ 

^ + 7^ + A2 H 2 kJ = -6^cmt4</'cm • (59b) 

In the regime (c) that we are considering now, the leading solution is obtained by neglecting 
the terms with time derivatives on the l.h.s. of (I59bp . This gives, 

where the expression for x is obtained using f l59ap . Clearly, these perturbations grow with 
time. As anticipated at the end of Sec. [3] this growth results from the low propagation speed 
of these modes. 

On top of the previous solution there is an oscillatory mode that corresponds to the 
solution of the homogeneous part of fl59b[) . 



1 f , ■ 
Uom cx - ■ exp ± z- 



t \ ^cmCeikct ^ 

However, it rapidly decays compared to (160|) and hence is irrelevant. 



Regime (d): subhorizon modes with k < yHakc. 

In this case the frequency of the modes is lower than the Hubble rate and we can neglect 
the third term on the l.h.s. of (I59bp . Then the perturbations are frozen and we obtain the 



same adiabatic solution flS^ as in the superhorizon case. It is important to stress though, 
that the modes we are discussing now have wave-lengths shorter than the Hubble size and 
are therefore observable. 



5.2 Corrections to the metric perturbations 

According to the iterative procedure described at the beginning of the current section, we 
shall substitute the expressions for the x and ^ perturbations obtained above into the Ein- 
stein's Eqs. f HOj) . This will give us the corrections to the observable quantities due to the 

^^An alternative derivation of this limit and Eqs. ([5^ involves diagonalization of the quadratic action (|26p . 



25 



presence of these fields. It is convenient to concentrate on Eq. fl40cp . Substituting in terms 
of and X from (140dl) we obtain 



2, a 

(2 + aB)M^ 



i) + mi) + (2-H + 

(61) 



2 + aB 

The r.h.s. represents the source for the correction t/''^^-' to the standard adiabatic behavior 
f l43|) of ip. Let us compare the two terms on the r.h.s. For superhorizon modes the second 
term dominates. Clearly, it gives a contribution of order 0(/3) to ip^^\ On the other hand, 
for subhorizon modes the first term becomes dominant. Indeed, it is enhanced compared to 
the second term by the ratio [k/T-Ly. Below we focus on the contribution of this term into 
ip^^^ and neglect the second term on the r.h.s. of ( 16T1) . 

Assuming as before that the relevant modes enter the horizon at the matter-dominated 
epoch, Eq. ( I6T1) is further simplified to 

This yields for -ip^-^^ the expression 

(i^0£',,v'«(,(n.ix(n). (62) 

where is the time of horizon crossing for the mode with given k. At this time, the adiabatic 
initial conditions imply tjj^^^t^:) ^ 0. 

Consider first the modes with kc < k. As it is clear from Fig. [H these modes spend all 
subhorizon evolution in the regime (6), where the field is given by ( |57|) . Inserting this into 
(!62l) gives, 

V'^^)~(/3 + A)-^^^-t^0e„.lnf . (63) 

X 

We see that for these short modes the correction to the gravitational potential remains para- 
metrically small, i.e. 0(/3 + A). Notice, however, that the numerical value of the logarithmic 
factor in (163|) can be quite large (of order 10), which can lead to a sizable effect if /3, A are 
not too small. 

Next we turn to the modes with ki/2 < k < k^, see Fig. [T] Due to fast growth of these 
modes in the region (c), the integrals in (162|) are saturated at the upper end of the integration 
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domain. Using (160!) we get an estimate for the value of ijj^^^ at presento: 



At /c ~ ki/2 we obtain 



where we have used Eqs. (1441) . fl45j) . We clearly see the parametric enhancement of the signal 
to the leve|3o(ya). 

Finally, we consider the modes with Hq < k < ki/2. Again, the integral ( 162|) is saturated 
at the upper end corresponding to the regime {d). Using the appropriate formula for x, 
Eq. ( !55l) . we obtain 

i^^'\t,) ^ + mt0?<Pcm ^ + A) (^-^^ <Pcm , 

where in passing to the last expression we have used Hq ^ 2/to- This expression has a 
maximum at the highest momenta of the interval, k ~ ki/2. Using ( 1451) one readily verifies 
that, up to a numerical factor of order one, this expression matches at ~ ki/2 with the 
estimate flMl) . 

To sum up, the analytical study of this section shows that the correction ip^^'^ to the 
present-day amplitude of the gravitational potential has a maximum of order y/a(f)cm at 
k ~ ki/2 with the falloff ip^^^ oc k'^, ip^^^ oc k~'^ on the two sides, and a logarithmic tale 
ip^'^^ ~ (/3 + A)0cm In k extending towards large momenta k > kc- This will be verified by the 
numerical results in the next section. 



6 Cosmological perturbations: numerical results 

In this section we report the results of numerical integration of Eqs. (l40l) . (14TI) . (l42l) . We 
consider the density fractions of radiation and cold matter today 

n^ = 5- 10'^ , n^rn = 0.25 . 

^^Strictly speaking, the formula (j62[) is applicable only during matter-domination. However, for the reasons 
discussed at the beginning of this section, we neglect possible corrections due to the start of the accelerated 
cosmological expansion. 

^^Recall that we assume the parameters a, (3, A to be of the same order. 
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Figure 2: Solution of the system (HOi) . (14T1) . (H2l) for two values of the momentum: k = 1.5Ho 
(upper panels), k = AHq (lower panels), where Hq is the present Hubble constant. Left panels 
show the dependence of the dark energy perturbations ^ on conformal time t. Right panels 
show the perturbations of the gravitational potential and the difference (0 — ip) between 
the two gravitational potentials. For comparison we also present the dependence of the 
gravitational potential in the standard ACDM cosmology. The parameters of the model are 
a = 0.02, P = 0.01, A = 0.01. For these values kc = IS.lSi^o, ki/2 = 3.89i/o, h/s = I-^GHq. 
Present time corresponds to to = S.SifJ'^ 



Correspondingly, the dark energy fraction is 



Qj^j^ = 0.75 
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Figure 3: The same as Fig. [2]for /c = 15ifo (upper panels), k = GOHq (lower panels). 



Let us emphasize again that for simplicity we neglect all effects related to baryon-photon 
coupling and merely include baryons into the cold matter fraction. Similarly, we neglect the 
neutrino masses. Still, within these simplifications we refer to the cosmological model with 
CC and without ^ and x perturbations as ACDM. Details of the numerical procedure are 
presented in the Appendix S Here we summarize the results. 

The time dependence of the fields the gravitational potential (j) and the difference 
(0 — ip) is shown in Figs. [2] and [3] for several values of the momentum k. The parameters of 
the model are taken to be 



a 



0.02, /3 = 0.01, A = 0.01 



C0 
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Figure 4: The ratio of density contrasts in the BCDM and ACDM models versus time for 
several values of the mode momentum k. The values of parameters are the same as for Fig. [2l 



This choice satisfies the PPN^I. BBN and gravitational radiation constraints [TTf 128] . The 
initial conditions correspond to the adiabatic mode, and, for illustration purposes, the initial 
value of (j) is normalized to 1 for every momentum. We also show for comparison the depen- 
dence of the Newton potential on time for the standard ACDM case (where ip = (p within 
our approximation). 

We see that perturbations of the gravitational potential are enhanced at late times com- 
pared to the ACDM case. The relative magnitude of enhancement is of order 10% which 
agrees with the analytic estimate ^/a 0.14 derived in the previous section. On the other 
hand, the difference between the two gravitational potentials (0 — tp) which is initially of 
order 10~^ rapidly decreases once the mode enters into the horizon, and as a consequence is 
at present negligible for k larger than ki/2. Note that the overall amplitude 10~^ for long 
wavelength modes agrees with the estimates {(p — ip) ^ a = 0.02. 

One expects that the enhancement in the gravitational perturbations will lead to the 
increase of structure growth rate. To illustrate this point we plot in Fig.Hlthe ratio S/6acdm, 
where 6 and Sacdm are the cold matter density contrasts in the case of the present model 

^^Recall that a — 2/3 case avoids the PPN bounds. 
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p. 
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/v 




"'1/2 yi-^ ivipL ) 


a 


2 • 10-2 


10-2 


10-2 


5.1 • 10-3 


1.3 - 10-3 


b 


2 • 10-3 


10-3 


10-3 


1.6-10-2 


2.3 - 10-3 


c 


2 • 10-^ 


10-^ 


10-^ 


5.0-10-2 


4.1 ■ 10-3 


d 


10-^ 





10-^ 


7.1 - 10-2 


4.9 - 10-3 



Table 1: Model parameters corresponding to the curves in Fig.|5]and Fig.[6]and the respective 
values of momenta kc, 

and ACDM respectively, 

S = (65) 

P[cm] 

As expected we observe the increase in the growth of structure at recent times. The relative 
effect is stronger for shorter modes and can be as large as 11% for our choice of parameters. 

The next plots show the comparison between the power spectra of perturbations in the 
0CDM and ACDM cosmologies evaluated at the present moment of time. Fig. O shows the 
relative differences of the power spectra for the gravitational potential and the cold matter 
density contrast 6, 

^'PW = p (Ta - 1' Mk) = p 77T - 1- 

These differences were computed for several values of the model parameters that are listed in 
Table [H We also present the corresponding values of the momenta kc and ki/2. Note that, as 
expected from the results of the previous section, the functions Aj(/c) peak around k ^ ki/2 
(different for the different sets of parameters) with a tail extending to larger momenta. 
In agreement with the analytical results, the value of the functions Aj at the peak scales 
approximately as ^/a while the tails are proportional to a. 

Finally, we plot in Fig. |6]the spectrum of the relative difference {(f) — ip) / (f) between the 
two gravitational potentials in BCDM. This difference is completely negligible at k larger 
than 0.01 Mpc-^. At smaller k it can be as large as a few percent for certain choices of 
parameters. 

We stress once again that all above results are valid in the case of longitudinal perturba- 
tions in the Einstein-aether model coupled to the 0-field. x this case corresponds to the 
longitudinal component of the aether perturbations. 
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0.001 
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Figure 5: Comparison between the power spectra for the gravitational potential (left panel) 
and the matter density contrast (right panel) in the GCDM and ACDM cosmologies. The 
curves correspond to the values of the model parameters listed in Table [TJ 




0.001 



k(hMpc'') 



0.01 



Figure 6: The difference between two scalar gravitational potentials for several choices of 
parameters listed in Table [H 



7 Conclusions 

In this paper we have proposed a dark energy model based on the idea that Lorentz invari- 
ance is not an exact symmetry of nature and is broken in the gravitational sector. As the 
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framework for Lorentz breaking we considered the theory of a unit time-hke vector field - 
aether. This can be either a general unit vector, in which case the setup corresponds to 
the Einstein-aether model, or it can be expressed in terms of the gradient of a scalar field 
- khronon - that defines a global time coordinate. Dark energy appears in this framework 
in a rather generic and simple way due to an additional massless scalar with an exact shift 
symmetry and arbitrary relevant couplings to the aether compatible with this invariance. An 
example of such scalar is the Goldstone boson of a global broken symmetry. The resulting 
model is a valid effective theory with a cutoff that can be only a few orders of magnitude 
below the Planck mass. Moreover, in the khronon case the model has a known candidate 
UV completion in the form of gravity with anisotropic scaling proposed by P. Hofava. 

We have shown that the model exhibits the property of self-acceleration. Namely, in the 
absence of the cosmological constant and any matter sources the model possesses two solu- 
tions corresponding to Minkowski and de Sitter space-times. The former solution is unstable 
and the presence of an arbitrarily small amount of matter destroys it. The cosmological 
evolution of a matter-filled universe is driven to the de Sitter attractor. The value of the 
effective cosmological constant in the de Sitter branch is determined by the lowest dimension 
coupling between the Goldstone field and the aether. Importantly, it is technically natural 
to assume this coupling to be small as it is protected from radiative corrections by a discrete 
symmetry. Thus if one were able to enforce somehow the vanishing of the vacuum energy and 
the existence of an (unstable) Minkowski vacuum, the value of the current cosmic accelera- 
tion would not present fine-tuning problems. However, we left the discussion of a mechanism 
that could lead to this cancellation of the vacuum energy outside the scope of this article. 

We analyzed the phenomenological consequences of the dark energy model proposed 
in this paper. To this aim, we supplemented the model with cold matter and radiation 
components. This particular realization was called OCDM. We showed that, barring fine- 
tuned initial conditions, the homogeneous cosmological evolution in the model is the same 
as in the universe with a cosmological constant. In other words, the effective dark energy 
equation of state is w = — 1. On the other hand, the non-trivial dynamics of dark energy 
reveals itself in the evolution of cosmological perturbations. We studied this evolution both 
analytically and numerically in the linear regime and found that the growth of perturbations 
is enhanced compared to the standard ACDM case. The enhancement is most prominent 
at very large scales of order a few gigaparsecs, but extends also to shorter scales. Another 
property of the model is the appearance of an effective anisotropic stress that leads at very 
large scales to the difference between the two gravitational potentials in the Newton gauge. In 
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principle, the previous effects can allow to discriminate OCDM from ACDM. Several groups 
have recently suggested a number of data analysis techniques to test the nature of dark energy 
and its possible deviations from the cosmological constant [361 EH [381 ES SQl El 1121 SSI Si] . 
It would be interesting to apply these techniques to the present case and work out the 
constraints on the parameters of the model following from the present data as well as their 
expected improvement by future experiments. 

Throughout the paper we have assumed a standard Lorentz invariant matter sector. For 
the Standard Model fields this follows from an overwhelming experimental evidence. How- 
ever, for the case of dark matter one could relax this assumption and allow it to interact 
directly with the Lorentz breaking fields. It would be interesting to study how these interac- 
tions may affect the evolution of cosmological perturbations. Potentially, this will allow to 
set bounds on the violation of Lorentz invariance in the dark matter sector. We leave this 
investigation for the future. 
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A The numerical procedure 

The complete set of equations for the cosmological perturbations is provided by the Ein- 
stein's equations ( HOll . equations (HTll for the khronon and Goldstone perturbations and the 
hydrodynamical equations fH2|) for matter (supplemented by the equation of state for each 
component). These equations are not independent. For the numerical solution we choose 
two out of the four Einstein's equations: Eqs. fl40cp and (140dl) . As indicated in the main 
text, we consider matter consisting of two components: radiation and cold matter. Only the 
former contributes into the pressure perturbation term in f l40cl) . This implies that to close 
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the system it suffices to consider the hydrodynamic equations (H2|) only for the radiation 
component. Introducing 

r _ ^PW 

Ph] 

these equations can be written as a single second order equation, 

Performing the Fourier decomposition, normalizing the present scale factor to one, a{tQ) = 1, 
and choosing the units such that the present Hubble parameter is equal to one, Hq = 1, we 
obtain the final system of ordinary differential equations to be solved numerically: 



^ + H{<j) + 2^) + {2H + H')<j) - 7r-—Bk\x + 2Hx) - tt^^h = , (66a) 

(j)-ij-/3{x + 2Hx) = , (66b) 
X + 2'Hx+ + -B)+ n^l + B) + a^cikl]x + aclkH 

= 4) + U{l + B)(t) + B^ , (66c) 

I + 2%^ + c|A;2f + ac%k'^x = -acl<p - 3aHc|0 , (66d) 
P 4P 

^1] + Y^h] + —<f> -4^ = 0. (66e) 

Here we have introduced the present radiation density fraction and defined ^ = ^//i^. 

A subtle point is the proper choice of initial conditions for the system (l66ll . These are 
fixed deep inside the radiation-domination stage when the modes are superhorizon. We 
consider the initial conditions corresponding to the adiabatic mode. The latter is regular at 
t —> 0. Thus we write for small t: 

= 0(o)+0(i)t, ^ = ^(0) + ^(1)^ ^ 5j^j=5(o) + 5Wt, (67a) 

x = x(°^t+x«^, e=e^°)^+e^^)^. (67b) 

z 2 6 

Expanding Eqs. ( 166|) at t — )■ we obtain the relations: 

= -20(°) , 6^ = ^^^^ , 
0(0) _ ^(0) _ 3^^(o) ^ Q ^ 0(1) _ ^(1) _ 2^^(i 



4y/a; 



x(°) = 0(°V2 , (3 + i?)x^'^ = (2 + 5)0^') + ^ 
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Here we have expanded the time dependence of the scale factor at the radiation domination 
epoch up to the subleading order, 

a = \ \l^t ^ 1 . 

Additionally, the initial data must satisfy the constraint following from the (00) Einstein's 
equatiorEf fl40ap . This gives, 

+ + ^ + + — 1 =0 



2 ^/n; V 2 

where 6 is the density contrast of cold matter, Eq. (!65|) . For the adiabatic mode the density 
contrasts are related: 

b = . 

Using this relation we obtain that all the coefficients in (l67|l are determined in terms of (p^'^^ 
Similar reasoning also gives the equations and initial conditions for ACDM. As we are 
interested in the ratios of various quantities in the GCDM and ACDM models we can choose 
arbitrary normalization for the only requirement being that this normalization is the 
same for the computations in both model^lj. In practice we choose (f)^^^ = 1. 
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